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We present results of a study of the decay J/^j) -^ uirj-K^-K^ using a sample of (225.2 ± 2.8) x 10® 
J/t/j events collected by the BESIII detector, and report the observation of a new process of J/^/j — > 
ajX(1870) in which X(1870) decays to Cq (980)7r^. The statistical significance of this process is 
larger than 7.2a. Signals for J/ip — > tj/i(1285) and J/ij) — ^ aj?7(1405) are also observed in rj-K^-K~ 
spectrum, with statistical significances much larger than 10a. 
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A new resonance, known as the X(1835), was first ob- 
served in the ry'Tr+Tr" mass spectrum of J/tp -^ 'yr]''K~^Tr~ 
by BESII [1[ and subsequently confirmed with a much 
higher signal significance by BESIII [2] . Several theoret- 
ical speculations have been proposed to interpret the na- 
ture of X(1835), including the pp bound state [3|-[a] that 
was first observed near the same mass in J/ip -^ jpp at 
BESII and confirmed by BESIII and CLEO 0, a sec- 
ond radial excitation of the r]' Q, and a pseudo-scalar 
glueball [9l-[ll|. In the lower mass region of the TyTr+Tr^ 
mass spectrum, around 1.4 GeV/c^, extensive studies 
[13 - [l4| have established the existence of the 77(1405), 
which has also been suggested as a candidate for a 
pseudo-scalar glueball |jj|. Experimentally, the study of 
the production mechanism of the X(1835) and 77(1405), 
e.g. searches for them in r]n~^Tr~ final states with other 
accompanying particles (w, 4>, etc.), are useful for clar- 
ifying their nature. In particular, the measurements of 
the production widths of these two states in hadronic 
decays of the J/ip and a comparison with correspond- 
ing measurements in J/ip radiative decay would provide 
important information about the glueball possibility [16| . 

In this letter, we present the results of a study of 
J/ip -^ uj-q-n^TT" . A structure around 1.8~1.9 GeV/c^ in 



the ?77r+7r^ mass spectrum is observed. This analy- 
sis is based on a sample of (225.2 ± 2.8) x 10^ J/ip 
event s |17| accumulated in the Beijing Spectrometer (BE- 
SIII) [1^ operating at the Beijing Electron-Position Col- 
lider (BEPCII) dl at the Beijing Institute of High En- 
ergy Physics. 

BEPCII is a double-ring e+e^ collider designed to pro- 
vide e+e" beams with a peak luminosity of 10^"^ cm~^s^^ 
at a beam current of 0.93 A. The cylindrical core of 
the BESIII detector consists of a helium-based main 
drift chamber (MDC), a plastic scintillator time-of- flight 
system (TOE), and a CsI(Tl) electromagnetic calorime- 
ter (EMC) that are all enclosed in a superconducting 
solenoidal magnet providing a 1.0 T magnetic field. The 
solenoid is supported by an octagonal flux-return yoke 
with resistive plate counter muon identifier modules in- 
terleaved with steel. The acceptance of charged particles 
and photons is 93% of 47r steradians, and the charged- 
particle momentum and photon energy resolutions at 
1 GeV are 0.5% and 2.5%, respectively. The BESIII de- 
tector is modeled with a Monte Carlo (MC) simulation 
based on GEANT4 [20l[2ll. 



Charged tracks in the BESIII detector are recon- 
structed using track-induced signals in the MDC. To op- 



timize the momentum measurement, we select tracks in 
the polar angle range | cos9\ < 0.93 and require that they 
pass within ±20 cm from the Interaction Point (IP) in 
the beam direction and within ±2 cm of the beam line in 
the plane perpendicular to the beam. Four tracks with 
net charge zero are required, and all tracks are assumed 
to be pious. 

Electromagnetic showers are reconstructed from clus- 
ters of energy deposits in the EMC. The energy deposited 
in nearby TOF counters is included to improve the re- 
construction efficiency and energy resolution. Show- 
ers identified as photon candidates must satisfy fidu- 
cial and shower-quality requirements, i.e. the showers 
in the barrel region (|cos0| < 0.8) must have a mini- 
mum energy of 25 MeV, while those from the endcaps 
(0.86 < I cos6i| < 0.92) must have at least 50 MeV. The 
showers in the angular range between the barrel and end- 
cap are poorly reconstructed and excluded from the anal- 
ysis. To suppress showers from charged particles, a pho- 
ton must be separated by at least 10° from the nearest 
charged track. The EMC cluster timing requirements 
are used to suppress electronic noise and energy deposits 
unrelated to the event. 

In the reconstruction of J/rp — > lotitt^tt^ , the uj is re- 
constructed via its TT+TT^TT*' mode, and the rj and ir'^ are 
reconstructed from 77 pairs. The vertex of all final state 
particles must be consistent with the measured beam in- 
teraction point. The sum of the four- momenta of all par- 
ticles is constrained to the known J/ijj mass and the ini- 



tial e 



+ r- 



three-momentum in the lab frame. The vertex 



and four-momentum kinematic fits are required to sat- 
isfy the quality requirements Xv ^ 1^0 and xlc < 5*-*' 
respectively. Further selections are based on the four- 
momenta from the kinematic fit. Photon pairs with an 
invariant mass satisfying AI^-^ G (524, 572) MeV/c^ or 
(122, 148) MeV/c^ are identified as rj or 7r° candidates. 
The rjir^Air combination with minimum xic i^ selected 
in the cases where more than one candidate satisfies 
the above requirements in an event. If there is more 
than one four-photon combination in the mass range of 
the 77 and tt", the assignment with the lowest value of 



XjjttO = \/-^v ~^ ^no ^^ used, where Pjj/^o are the pulls de- 
'^/^o — — ^^^ — 2/z_^ Here a^^ is the rj/n'^ mass 



fined as P, 

resolution determined from data. 

After the application of the above requirements, the 
scatter plot of M^(^^) versus M(^(^+^^^o) (shown in 
Fig. [1]) shows a clear cluster in the J/ip -^ ujr]Tr~^TT~ sig- 
nal region denoted by the rectangle in the center of the 
plot. To determine which tt+tt^ pair originates from the 
cj, \M^+„-^o —rrit^l is minimized among the possible com- 
binations of the selected charged pious, and required to 
be less than 28 MeV/c^. 

With all the selection criteria applied, the mass spec- 
trum of rin~^'!T~ is shown in Fig. [IJa). In the lower 
mass range, in addition to the well-known ry' peak, two 
other structures are observed; these are inferred to be 
the /i(1285) and 77(1405) based on the fit results dis- 
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FIG. 1: Scatter plot of M, 



v(-ii) 



versus M 



(ir-+ 



-ttO)- 



The 



rectangle in the middle shows the signal region defined as 



\M^ 



< 28 MeV/c^ and \Mj.y -m^| < 24 MeV/c^ 



cussed below. There is an additional structure located 
around 1.87 GeV/c^ that we denote as X(1870). The 
r/TT^ mass spectrum for these events shown in Fig. [2I^b), 
reveals a strong ao(980) signal. The ryTr+Tr" mass spec- 
trum for events where either M(r]7T'^) or M(r]Tr~) is in a 
100 MeV/c^ mass window centered on the ao(980) mass is 
shown in Fig.l^Jc). The 777r"'"7r~ mass spectrum for events 
with both Af (777r+) and Af (j/tt^) outside the ao(980) sig- 
nal region is shown in Fig. [2][d). A comparison between 
of Fig. ^c) and Fig. HJd) indicates that the /i(1285), 
77(1405) and A:(1870) all primarily decay via the ao(980)7r 
channel. 

To ensure that the observed /i(1285), 77(1405) and the 
structure around 1.87 GeV/c^ originate from the pro- 



cess of Jf-ip 



TT^ rather than peaking back- 



grounds, potential background channels are studied us- 
ing both data and MC samples. The non-w and/or non- 
ao (980) processes are estimated by the events in the two- 
dimensional mass-sidebands of w and ao(980) illustrated 
by the dashed and dotted boxes in Fig. [31 where the 
difference of the background shape in the signal region 
and side bands due to the varying phase space is taken 
into account by multiplying the rjiT^TT' mass distribu- 
tion by a correction curve. In practice, this contribution 
is evaluated by the weighted-sums of horizontal and ver- 
tical side bands, with the entries in the diagonal side 
bands subtracted to compensate for the double count- 
ing of background components. The weighting factors 
for the events in the horizontal, vertical and the diagonal 
side bands are measured to be 0.48, 1.58 and 0.76, respec- 
tively, which are determined from the results of a two- 
dimensional fit to the mass spectrum of M^{7r^TT^Tr'^) 
versus Mag[gso)(jl''^)- Here the two-dimensional Proba- 
bility Density Functions (PDFs) for J/ip — >■ a;ao(980)7r, 
UJ but non-ao(980), non-w but ao(980), non-cj and non- 
ao (980) processes are constructed by the product of one- 
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FIG. 2: Invariant mass distributions for the selected candi- 
date events: (a) and (b) are the invariant mass spectra of 
rj-K^Tv' and jjtt after the apphcation of all the event selec- 
tion criteria; (c) is the r)-K'^-K~ mass spectrum for events with 
an ao(980) in the -q-K^ final state; (d) is the tj-k'^'k' invariant 
mass distribution for events with no ao(980) in the -q-K sys- 
tem. The histograms in (a) and (c) are the phase-space MC 
events of J /ip — >■ Lori-R^-K~ after the same event selection and 
with arbitrary normalization. 



dimensional functions, where the resonant peaks are 
parametrized by Breit-Wigner functions and the non- 
resonant parts are described by floating polynomials. A 
MC sample of two million events of the phase space pro- 
cess of J/V' — > 7r+7r~7r''7y7r"^7r~ is used to determine the 
correction curve to account for the varying phase space. 
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The background channel J/i) -^ 5i(1235)ao(980), 
where the 6i(1235) decays to ujtt and ao(980) decays to 
rjTT, is studied by performing a two-dimensional fit to 
the M(u!tt) versus Af (t^tt) mass distribution with two- 
dimensional PDFs defined in similar fashion. We also 
studied an inclusive MC sample of 2 x 10^ J/ip decays 
generated according to the Particle Data Group (PDG) 
decay table and Lund-charm model ^^. The reliability 
of the background estimation method described above 
is validated, and no background-induced peaks are ob- 
served around 1.87 GeV/c^. 

Figure S] shows the results of a fit to the ryTr+Tr^ 
mass spectrum where either rjn^ or rjiT^ are in the 
ao(980) mass window. Here the three signal peaks are 
parametrized by efficiency-corrected Breit-Wigner func- 
tions convolved with a Gaussian resolution function, the 
width of which is determined from signal MC samples 
and fixed in the fit. The background consists of three 
components, namely contributions from non-w and/or 
non-ao(980) processes, J/ip -^ 6i(1235)ao(980) events, 
and non-resonant a;ao(980)7r processes. The background 
shapes and numbers of events for the non-o; and/or non- 
ao(980) processes are determined from the events in the 
two-dimensional side bands as discussed above, and fixed 
in the fit. For the J/ip -^ 5i(1235)ao(980) component of 
background, the background shape is fixed to that of the 
phase space MC samples whereas the number of events 
is extracted and fixed to the result of a two dimensional 
fit to the WTT versus rjir mass distributions. The contri- 
bution of the remaining non-resonant wao(980)7r process 
is described by a smooth floating polynomial function. 
The mass, width, and the product branching fractions 
obtained from the fit are summarized in Table [H For the 
/i(1285) and 77(1405), the measured mass and width are 
in agreement with PDG values [23| . A conservative esti- 
mate of the statistical significance of the X(1870) signal 
is determined by the lower limit in the change of — 2 ln>C 
obtained from the fits with and without the assumption 
of a X(1870). With aU the factors in the fit varied, the 
smallest change in the — 21n£ is 60.1, corresponding to 
a significance of 7.2a. The same procedure is applied to 
the /i(1285) and 77(1405) signals, and the significances 
are determined to be much higher than 10 a. 



TABLE I: Summary of measurements of the mass, width 
and the product branching fraction of B{J/il) — >■ ljX) x 
B{X -^ a±(980)7r'F) x B(aJ(980) -^ t/tt*) where X repre- 
sents /i(1285), »7(1405) and X(1870). Here the first errors 
are statistical and the second ones are systematic. 



Resonance 



/i(1285) 
Jj(1405) 
J>C(1870) 



Mass (MeV/c^) Width (MeV/c^) B (lO^*) 



1285.1 ± l.O! 



+ 1.6 
-0.3 

1399.8 ± 2. 2lo- J 
1877.3 ± 6. 3I7-4 



22.0 ±3.1 

52.8 ±7.6 

57 ± I2I 



+2.0 
-1.5 
+0.1 



i.25±o.ior^;^^ 

1.89 ±0.21+°;^! 
1.50±0.26lo;36 



FIG. 3: Definition of the signal and two-dimensional side 
bands. 



The systematic errors on the measurement of the mass 
and width parameters are primarily due to the uncer- 




^777r"i"7r 



FIG. 4; Results of fits to the M{riTT'^TT~) mass distribution 
for events witli either the 7771 "'" or r]n~ in the ao(980) mass 
window. The dotted curve shows the contribution of non-o; 
and/or non-ao(980) background, the dashed line also includes 
the contribution from J/ip -^ 6i(1235)ao(980), and the dot- 
dashed curve indicates the total background with the non- 
resonant J/ip — >■ ojoq (980)7r^ included. ^ jd.o.j is 1.27 for 
this fit. 



tainty in the mass spectrum fitting. In detail, the fit 
range, background estimation method, number of back- 
ground events, and the background parametrization are 
varied to decide the uncertainty from the background es- 
timation and fitting as a whole. We also include the sys- 
tematic errors determined from the input/output checks 
based on the analysis of full-reconstructed MC sam- 
ples, in which the input parameters are set according 
to the final results and the background is represented 
by the background channels seen in the inclusive MC 
sample. For systematic errors originating from the po- 
tential structure around 2.2 GeV/c^ and the multiple- 
event candidate selection, we re-fit the mass spectrum 
of 777r+7r~ with the inclusion of a X(2120) resonance as 



recently reported by BESIII [2| in the decay channel of 
J/^/i — >■ 77y'7r+7r~, and all the valid multiple-entry can- 
didates kept in order to estimate the uncertainty due to 
these two sources, respectively. With the numbers from 
all the sources above combined quadratically, the sys- 



tematic errors on the mass and width parameters are 
determined as shown in Table |ll 

The systematic errors on the branching fraction mea- 
surements are also subject to errors of the number of J/V' 
events ^^] , the intermediate branching fractions [2J| , the 
data-MC difference in the tt tracking efficiency, the pho- 
ton detection efficiency, the kinematic fit, the signal se- 
lection efficiency of jy/Tr", the simulation of the line shape 
of ao(980) |2j], and the angular distributions due to 
different possible spin-parity hypotheses. Combined in 
quadrature with the influence from the mass spectrum 
fitting, the systematic errors on the product branching 
fraction for the /i(1285), r/(1405) and X(1870) are sum- 
marized in Table m 

In summary, we present a study of the J/V" — > w?77r+7r~ 
decay channel and report the first observation of a new 
process of J/t/i -h- ^^(ISTO) in which X(1870) decays 
to Aq (980)77^, with the signal significance estimated to 
be 7.2cr. In the lower mass region of t^tt+tt" mass spec- 
trum, the /i(1285) and ?7(1405) are also clearly observed 
with statistical significances much larger than lOcr. The 
measurements of the mass, width, and product branch- 
ing fraction of B{J jii -^ cjX) x B{X -> a^{980)n^) x 
S(aQ (980) -^ ?/7r*) for the three resonant structures are 
summarized in Table [H wherein the branching fractions 
for the /i(1285) and 7/(1405) are measured for the first 
time. Whether the resonant structure of X(1870) is due 
to the X(1835), the 772(1870), an interference of both, or 
a new resonance still needs further study such as a par- 
tial wave analysis that will be possible with the larger 
J/ 7/; data sample that is anticipated in future runs of the 
BESIII experiment. For the 77(1405), the product branch- 
ing fraction of its hadronic production is measured to be 
smaller than that for its production in the radiative J/tp 
decays |23j |. 
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